Squeezing and entanglement of atomic beams 
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We propose and analyze a scheme for generating entangled 
atomic beams out of a Bose-Einstein condensate using spin- 
exchanging collisions. In particular, we show how to create 
both atomic squeezed states and entangled states of pairs of 
atoms. 

PACS numbers: 03.75.Fi, 03.67.Hk, 42.50.-p, 03.65. Bz 
The recent experimental achievement of Bose-Einstein 
condensates has raised a lot of interest since it may 
lead to important applications j^-^ . Some of these ap- 
plications are based on the analogies between an atomic 
condensate and a single mode optical field. For example, 
atoms in the condensate can be outcoupled to produce a 
coherent atomic beam Q similar to a laser beam. In this 
letter we build on these analogies to show how to pro- 
duce entanglement between atoms in different internal 
states similar to the one created for photons with dif- 
ferent polarizations by parametric down conversion [^,0 . 
In particular, we analyze a physical situation which gives 
rise to: (i) a beam of atoms in a broad-band two-mode 
squeezed state (a continuous variable entangled state) 
with respect to two internal levels; (ii) a pair of outgoing 
atoms in an effective maximally entangled state in a two- 
dimensional Hilbert space. The physical mechanism re- 
sponsible for this processes is spin exchanging collisions, 
where two atoms of the condensate interact to create two 
correlated atoms in two different internal states. 

We consider a Bose-Einstein condensate confined in an 
optical trap and in some internal level |0). Two atoms 
in the condensate can collide to create a pair of atoms in 
two other internal levels 



2|0)^| + 1) + |-1). 



(1) 



The situation we have in mind is illustrated in Fig. I. 
Levels |0, ±1) could correspond to the hyperfine Zeeman 
levels |F = 1, M = 0, ±1) of an Alkah atom. In that 
case, the selection rules would prevent other collisional 
processes to occur. We assume that the process (0) can 
be switched on and off by changing some external param- 
eter. For example, the condensate level could be shifted 
in a time-dependent way by an external field so that en- 
ergy conservation effectively allows or inhibits the pro- 
cess in Eq. (|l|). This level shift could be accomplished, 
for example, by an off-resonant microwave or laser field 
with an appropriate polarization. We also assume a one- 
dimensional situation where the trapping potential along 
the transverse direction is sufficiently tight so that the 
atomic motion is frozen along the y, z directions. We will 
choose different trapping potentials along the x direction 



in order to illustrate our ideas of how to create squeezed 
atomic states and entangled pairs. In both cases, the 
potential will be identical for the atomic levels | ± 1) 
and such that the atoms in those levels can escape from 
the trap, in order to facilitate measurements with them 
without being affected by the atoms in the condensate 
level |0). For instance, this could be obtained if the op- 
tical trap is made by a strong laser along the x direction 
crossed by a weak laser along y. In this configuration, 
if the energy shift between the |0) and | ± 1} levels is 
larger than the trap depth induced by the strong laser 
and smaller than that by the weak laser, the \Mp — ±1) 
atoms will be free to move along x but will be bound in 
the other two directions. 
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FIG. 1. Atomic configuration: The condensate is in the 
Hyperfine Zeeman level |_F = 1, Af = 0), which is shifted with 
respect to the \F = 1, M = ±1) states. 

This situation is described by the following second 
quantized Hamiltonian [h ~ \) 



2m 



(j)i{x)dx 



+ r g{x,t) Ul,{x)^l,{x)e-^^^''+h.c) dx, (2) 

where, (j)±i are quantized fields describing the atoms in 
the internal levels | ± 1) satisfying the commutation re- 
lation [(^i(a;, i), 0](a;', i)] — 5ij5{x — x'). In Eq. (2), the 
first term describes the kinetic and potential energy of 
the atoms in the internal levels | ± 1). The second term 
describes the creation (or annihilation) of two atoms in 
those levels due to a collision between two atoms in the 
condensate. The condensate is described by a macro- 
scopic wave function and /i is the corresponding chemi- 
cal potential. The function g{x,t) is proportional to the 
s-wave scattering length and the condensate density. We 
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have included the time-dependence explicitly to account 
for the depletion of the condensate, as well as to take into 
account the change of the external parameters which al- 
low to control the process (|l|). On the other hand, in 
writing Eq. (2) we have ignored the collisions among the 
atoms in levels | ± 1} as well as the quantum fluctua- 
tions of the condensate, which is valid if the number of 
atoms in the condensate is large and in the other levels is 
small 1^. There are other coUisional terms, like ^\4>i(t>Q(t>Q 
{i = ±1), which are allowed by the coUisional selection 
rules and are of the same order in the number of con- 
densate atoms as the second term in the Hamiltonian. 
However, they can be included as an effective potential 
for the atomic beam fields. In the following, we under- 
stand V (x) as a renormalized potential which includes 
such coUisional terms. 

The quadratic Hamiltonian (2) gives rise to the follow- 
ing linear Heisenberg equations 

tdt(t)±iix,t) = --^(l)±i{x,t) +V{x) <l>±i{x,t) 
Am 

+g(x,t)0^i(a;,i)e-^2^*. (3) 

In the Heisenberg picture, the initial state of the atoms 
1 4*) in levels | ± 1) is the vacuum, since we assume that 
initially (at time t — + — cxd) all the atoms are in the con- 
densate, that is, (t)i,2{x, — oo)|^') = 0. Thus, the coupled 
equations describe the generation of atoms in levels 
± 1) out of the vacuum. Given a potential V{x) and 
the function g{x,t) one can solve them numerically in 
the same way as one solves the Bogoliubov-de Gennes 
equations for the time-dependent excitations of a con- 
densate 1^. Thus, the solutions can be written in terms of 
time-dependent Bogoliubov transformations, which im- 
plies that there will be correlations between the atoms 
created in levels | ± 1). This is precisely the physical 
origin of the entanglement that will be described in the 
following. 

We will now study two limiting situations: (i) The 
Squeezing limit, in which the typical time for the atoms 
to leave the trap is much smaller than the one related to 
the coUisional process, so that many atoms accumulate 
in the levels | ± 1) before they escape from the trap, (ii) 
The qubit Entanglement limit, where the atoms leave the 
trap before a new pair of atoms is created. The first sit- 
uation is analogous to the squeezed light generation by 
non-degenerate parametric down conversion Q , whereas 
the second case is similar to the one in which pairs of 
polarization entangled photons are created 

Let us first study the situation of squeezed atomic 
beams. We consider a simple model for which we can 
obtain the properties of the atomic correlations analyti- 
cally. The conclusions drawn from this model are qual- 
itatively valid for more complicated situations, in which 
one has to rely on numerical calculations. We take a 
potential V{x) = for a; > and infinite otherwise, so 



that 0±i(O, t) = 0, and assume that the condensate wave 
function is such that g{x,t) = go{t) for < x < a and 
zero otherwise. The function go{t) is switched on and 
off in a time of the order 7"^. This time can be re- 
lated, for example, to the typical depletion time of the 
condensate. We will consider separately the two spatial 
regions: (I) x > a] (II) < a: < a and then connect them 
via the requirement that the field operators and the first 
derivatives have to be continuous. 

For X > a (I) we can write the solutions of Eq. (||) as 

(^±i(x,t)= J (i±i(cc;)e-^'=('")"+B±i(w)e''=(")")e-'-*dLj, 

(4) 

where fc(aj) = ^^miolh. Scattering theory assigns 
the operators Aj^-y{lo) and B±i{lo) a definite physi- 
cal meaning. They are annihilation operators of par- 
ticles in incoming and outgoing plane wave modes 
with velocity =F\/2?iw/m, respectively. The condition 
(/>±i(x, — oo)|\['} = is then translated into A±i(w)|5') = 
0. The physical interpretation is that initially there are 
no input atomic beams, so the input modes should be 
in the vacuum state. The output modes -B±i(w), deter- 
mined by the inputs A±\ [lo) and the dynamics in the con- 
densate region, are directly related to measurable quanti- 
ties. The state of the output components can be detected 
by velocity-selective light imaging In order to deter- 
mine the state of the output modes for vacuum inputs, 
we need to solve the dynamical equation (|^) in the con- 
densate region < x < a. 

For this purpose, it is convenient to take a Fourier 
transformation of Eqs. (^, obtaining a coupled set of 
equations for "I>±i(x, A) defined through 

/oc 
$±i(a:,A)eT'^*dA. (5) 
-00 

Due to the fact that go (0 is time dependent, the Fourier 
components <I'+i(x, A) are correlated with $_i(x. A — 
for a range of frequencies lo. For applications, however, it 
is desirable to have pure entanglement between two mea- 
surable modes, which in our case are the output Fourier 
components. In fact, that can be obtained in the limit 
J go, where go is the maximum value reached by 
go{t). As shown below, the bandwidth of <I>±i(x, A — ui) 
is roughly determined by the coupling rate go, which is 
much larger than the width of the Fourier transform of 
go{t). This means that the Fourier transform of go(t) can 
be replaced by gQd{u!). We call this approximation the 
steady output condition, since it corresponds the phys- 
ical condition that the atomic loss in the condensate is 
negligible before we get a steady output. Imposing the 
boundary conditions allows us to express the outgoing 
operators in terms of the ingoing ones as a Bogoliubov 
transformation 
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B±i(a* ± A) = «±i(A)i±i(^ ± A) + p±i{A)Al,{fi T A), 

(6) 

where the coefficients a, P can be determined by solv- 
ing the corresponding scattering equations. Note that to 
keep the commutation relations these coefficients satisfy 
the general requirements |q;±i(A|^ — |/3±i(A)p = 1 and 
a+i(A)/3_i(A) - a_i(A)^+i(A) = 0. From Eq. (|), we 
see that the outgoing modes -B±i(/i ± A) are in a pure 
two- mode squeezed state with the squeezing param- 
eter TA given by 



tanh (rA) 



1/3+1 (A) I _ |/3_i(A)| 
|«+i(A)| |a_i(A)|- 



(7) 



The dependence of the squeezing r/\ on the detuning A 
determines the squeezing spectrum. Note that a pure 
two-mode squeezed state is an ideal continuous variable 
entangled state, with the entanglement characterized by 
the squeezing parameter [|Io|. Continuous variable en- 
tangled states have many application in recent quantum 
information protocols . 

To simplify the expression for the squeezing parameter 
ta, we assume fJ.^ go- This can be achieved in practice 
since /j, is adjustable by changing the shift of level |0). In 
this case, the squeezing parameter can be written in the 
following simple form 



TA = jarctanh {tanh(26') sin [(fc+ — fc_)a]} 



(8) 



where 9 = arctanh[[(A/go)^ + 1]^^^ - A/.go], and k± = 
[2mn/h ± (2m.go/n)[(A/.go)^ + l]^/2]i/2_ ' xhe solution 
reveals some interesting properties of this interac- 
tion. First, let us consider a vanishing detuning A = 0, 
that is, we look at the squeezing tq between the output 
modes and _B_i(/i). The parameter vq is given by 

ro = jarctanh [sin (go^)] |i where t — 2a/ y^2hfi/m is the 
transmission time of the input atomic beam with velocity 
^J2?iiji/m in the region < x < a. If the dimensionless 
interaction coefhcient k = 170^ — 7r/2, we have infinite 
squeezing and infinite output atomic fiux. This simply 
means that the approximation of negligible atomic loss 
for the condensate has broken down at this point. So, 
similar to the non-degenerate parametric down conver- 
sion in the optical case our system has a working 
threshold given by k = 7r/2 -I- nir. The system should op- 
erate not very close to the threshold to get steady output 
of entangled atomic beams. 

Next, let us look at the squeezing spectrum. The 
squeezing ta versus the dimensionless detuning A/(7o 
and the interaction coefficient n = got is shown in Fig. 
2. From the figure, we see that we have a broadband 
two-mode squeezed state with the squeezing bandwidth 
roughly determined by g^. The steady output condition 
requires go ^ 7- In our case, the atomic loss is mainly 
caused by the output coupling. From Eqs. (0) and (|^), 
the loss rate can be estimated as 7 ~ 2f;osinh (ro) /Nq, 



where Nq is the total atom number in the condensate. 
Even with a high peak squeezing ro, the steady output 
condition can still be easily attained. It is also interest- 
ing to note that one can control the transmission time 
t by changing the level shift to get a large peak squeez- 
ing ro, and we know that the steady output condition 
does not put a stringent requirement on the obtainable 
value of ro. Thus, in this system in principle one can 
get a much larger squeezing and therefore a much larger 
entanglement than in the optical system. As an exam- 
ple, a conservative estimate gives g^ ~ 20kIIz, a ^ 3/im, 
iJ — y^2hfi/m ^ 9cm/s, corresponding an output flux 
about 680 atoms/ms, we have a very large squeezing 
ro ^ 2, which is not yet obtainable in current optical 
systems. The advantage of large obtainable squeezing 
in this system is due to the fact that we have a strong 
nonlinear interaction caused by the collisions with the 
condensate. 

FIG. 2. Squeezing parameter r versus dimensionless detun- 
ing A /go and interaction coefficient k = got 

We emphasize that despite the simple model for the 
potential and condensate shape, we expect that all these 
features will be present for more realistic models. In fact, 
if the potential is asymptotically flat we can write Eq. (Q) 
in that region, so that under the steady output condition 
we will obtain Eq. (^) with different coefficients a and 
/3. For demonstration of these features, we suggest a 
three-step experiment. First, one can demonstrate the 
existence of a threshold by controlling the velocity of the 
output beams. For a certain velocity of the beams, the 
spin relaxation rate of the condensate will increase dra- 
matically, and that gives the threshold value. Secondly, 
one can measure the squeezing spectrum by the veloc- 
ity selective light imaging Finally, one can demon- 
strate the entanglement between the atomic beams by 
atomic homodyne detection This can be achieved 

by a atomic beam splitter and a local oscillator provided 
by an atom laser, which can be outcoupled from the same 
condensate. 

Now, we will study the situation in which pairs of 
atoms are created sequentially. We will show that, as 
in the case of photons if we post-select the mea- 

surement results, the corresponding internal state of the 
atoms is effectively maximally entangled. For this pur- 
pose, we assume the condensate is located at the region 
—a < X < a, and the potential V (x) is symmetric and 
independent of the internal states so that the | ± 1) atoms 
have the same probability of going to the regions x < —a 
and X > a. Two atomic detectors are placed, one on 
the left {x < —a) and one on the right {x > a) of the 
condensate. The coupling g{x,t) in the Hamiltonian (j^) 
is assumed to be sufficiently small such that there are 
at most one pair of atoms generated in each detection 
interval. Using perturbation theory (in the Schrodinger 
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picture) we can obtain from the Hamiltonian ([ 
fective atomic state for each detection interval: 



the ef- 



^(t)) = / fix,y,t)4>l^ix)4>l,{y)dxdy\vac), (9) 



where f{x, y, t) can be easily calculated and we have ne- 
glected the vacuum component since it has no influence 
on the measurement results. After a time the atomic 
pair leaks out of the condensate, and the wave func- 
tion f{x,y) = f{x,y,to) ^ ioT —a < x,y < a. We 
can decompose the wave function in the form / {x,y) = 
Jlr, [x, y) -f fRL {x, y) + Jll [x, y) + Jfir [x, y), where 
fLR{x,y) is defined to be equal to f {x,y) if a; < —a 
and y > a, and to be zero elsewhere. Other components 
are defined in a similar way. So the state \^ (to)) is de- 
composed into four components, with definite physical 
meaning for each component. For instance, the compo- 
nent /ll (x, y) represents both of the atoms go to the left 
side. Now we project the state onto the subspace where 
there is one atom at each side. This projection can be eas- 
ily achieved in experiments by post-selections of the mea- 
surement results, similar to many optical experiments in- 
volving spontaneous parametric down conversion |0,|l2|. 
After the projection, we only have two components in 
the effective state l^'e//) (the state selected by the mea- 
surement). The potential is independent of the internal 
states of the atoms, so we have /lr {x,y) — Jrl iy,x). 
With this condition, the effective state has the form (not 
normalized) 



|*e//) - J fLR{x,y) (10) 

(i>+i (x) 4>li (y) + (x) 4>l^ (y) dxdy\va.c), 



LR 



LR 



which can also be written as — 1) 
with a different notation. Thus with post-selection of 
the measurement results, we get an effective maximally 
entangled state between the atomic pair, and this state 
should have many applications in the field of quantum in- 
formation, such as measurement of Bell inequalities with 
massive particles fl^ . 

In summary, we have analyzed a scheme for generat- 
ing both entangled atomic beams and entangled atomic 
pairs. We have shown that we can get pure continuous 
entangled state with a large entanglement and effective 
maximal qubit entanglement with post-selection. The 
generated pure atomic entanglement can be directly used 
in the demonstration of many interesting quantum infor- 
mation protocols p4[ . 
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Note: Upon completion of this work, we have become 
aware of the very recent preprint quant- ph/0007012| by 
H. Pu and P. Meystre, which discusses similar ideas on 
creating squeezed atomic states. 
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